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The elastic energy functional of a system of discrete dislocation lines is well known from dislocation 
theory. In this paper we demonstrate how the discrete functional can be used to systematically derive 
approximations which express the elastic energy in terms of dislocation density-like variables which 
average over the discrete dislocation configurations and represent the dislocation system on scales 
above the spacing of the individual dislocation lines. We study the simple case of two-dimensional 
systems of straight dislocation lines before we proceed to derive energy functionals for systems of 
three-dimensionally curved dislocation lines pertaining to a single, as well as to multiple slip systems. 

We then illustrate several applications of the theory including Debye screening of dislocations in 
two and three dimensions, and the derivation of back stress and friction stress terms entering the 
stress balance from the free energy functionals. 

PACS numbers: 46.50.-|-a, 62.20.M-, 62.20.mm, 64.60.av 


INTRODUCTION 


Any dislocation density based theory of dislocation dynamics under stress, and thus of plasticity, must of necessity 
consist of two parts. The first consists in a density based description of dislocation kinematics, i.e., of the way 
dislocations as curved and connected lines move in space. The second consists in a density based description of 
dislocation energetics which allows to derive, via the powerful mathematical tools of variational calculus, the driving 
forces for dislocation motion. If one wants to go beyond the standard concepts of linear irreversible thermodynamics, 
a closed theory may require a third ingredient which provides the, in general non-linear, connection between driving 
forces and dislocation fluxes. 

As to describing the kinematics of dislocations on scales above the dislocation spacing where individual dislocation 
lines can no longer be resolved, significant progress has been made in recent years. In particular, the important 
question how to correctly describe the coupled kinematics of ’statistically stored’ and geometrically necessary dislo¬ 
cations for three-dimensional (3D) dislocation systems has been addressed in several works [H-Q- The progress in 
dislocation kinematics calls for a matching effort to develop averaged, density based descriptions of the energetics of 
dislocation systems. Such an effort needs to consider both excess dislocations associated with the spatial average of 
the classical dislocation density tensor (often termed geometrically necessary dislocations) and so-called statistically 
stored dislocations of zero net Burgers vector which dominate plasticity in the early stages of deformation. 

For dislocation systems described as systems of discrete lines, expressions for the associated elastic energy functional 
and interaction stresses have been provided in the classical works of dislocation theory, see e.g Q and, for overwiew, 
the textbook of Hirth and Lothe [9j. Recent progress has focused on methods to regularize the elastic singularity 
and associated diverging energy density in the dislocation core, using either nonlocal elasticity theories (e.g. [l^ or 
continuous Burgers vector distributions [ll|. Starting from the energy functional of a discrete dislocation system, one 
may then develop an energy functional for dislocation densities through appropriate averaging procedures. This is 
the approach pursued in the present investigation. 

There are alternative approaches. In the traditional spirit of constitutive modelling in continuum mechanics, 
numerous authors have introduced free energy functionals which depend on dislocation density-like variables where 
the related functional dependencies were assumed in an ad-hoc manner, see [l^ as one example of many. Usually, 
little attention is paid to the question how, if at all, these functionals can be derived from properties of the underlying 
discrete dislocation systems. We hnd this method of reasoning of little relevance to our own work, however, we note 
that our results can serve as a benchmark to assess whether or not the structure of various free energy functionals 
used in the literature is adequate for representing actual dislocation systems. Groma and co-workers derived equations 
for dislocation density evolution in 2D by direct averaging of the discrete dynamics l3-15| and then proceeded to 
design, via educated guess, free energy functionals which are consistent with this averaged dynamics [l^, or which are 
directly designed to reproduce average behavior in discrete dislocation dynamics simulations [ 13 . In the present work 
we pursue a modified but related approach - instead of systematically averaging the dynamics, and then matching it 
with free energy functionals obtained via educated guess, we apply systematic averaging to directly obtain an energy 
function. 
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Recently, several authors have attempted to derive free energy functionals for dislocation systems using thermody¬ 
namic formalisms, 


see e.g. 


18- 


In the opinion of the present author, thermodynamic approaches need to deal 
with the fundamental problem that dislocations do not exist in thermal equilibrium, and some authors may not be 
sufficiently aware of the implications of this basic fact. To illustrate the problem, let us compute the energy associated 
with a single Burgers vector of dislocation line length. This energy is of the order of fj,b^ where b is the length of the 
Burgers vector and /r is the shear modulus of the material. For typical materials parameters of copper, this energy 
amounts to about 4 eV. To provide matching thermal energies one would need to consider temperatures of the order 
of 40000 K, which is more than one order of magnitude above the melting temperature. Kooiman et. al. [l^ make 
the same observation when they notice that the coupling constant which gives the relative magnitude of elastic to 
thermal energies is, for dislocation systems at room temperature, of the order of 100. Hence, thermal effects and thus 
classical thermodynamic entropy are practically irrelevant to dislocation systems, whose driving forces derive almost 
exclusively from the internal (elastic) energy. Attempts to evaluate properties such as the range of correlations in 
dislocation systems using thermodynamic formalisms have produced interesting results, see e.g. the remarkable work 
of Limkumnerd and Van der Giessen who use a Langevin-type approach to evaluate correlation functions, with re¬ 
sults that are consistent with the results of discrete dislocation dynamics (DDD) simulations. However, these authors 
recognize that the fluctuation magnitude that must be assumed to achieve such agreement is orders of magnitude 
above the level of thermal fluctuations even at the melting temperature, hence it represents some kind of ’’effective 
temperature” - in fact a fit parameter that needs to be adjusted to make the range of dislocation correlations equal 
to a few dislocation spacings as observed in simulations. An effective temperature is also introduced by Groma et. 
al. [l^ in their free energy expression, and these authors make the same observation. These approaches may be 
considered implementations of the general suggestion by Berdichevsky to consider microstructural disorder in terms 
of an effective microstructural entropy and associated temperature [22j| . 

The bottomline is that, in order to explnin why for instance the screening radius in a dislocation system is of the 


order of several dislocation spacings [1^ l23| , one needs to introduce effective temperatures that cannot be related to 
standard temperature and thus, effectively, constitute phenomenological fit parameters. 

So why do we find, in plastically deformed crystals, densities of ’statistically stored’ dislocations (dislocations 
which have zero net Burgers vector and might thus annihilate) which are as high as 23|, when equilibrium 

thermodynamics requires this density to be zero? The answer is simply that statistically stored dislocations exist 
because kinematic constraints prevent them from annihilating. One can understand the problem best by considering 
dipoles consisting of dislocations of opposite sign moving on parallel slip planes: As long as the interaction is not 
sufficient to overcome the high energy barrier that prevents dislocations from leaving their slip planes, the dislocations 
will form a dipole with a width that is dictated by the slip plane spacing. Furthermore, dislocations are likely to 
become trapped in the first local energy minimum close to their initial position, and thermal energies (which are 
orders of magnitude less than dislocation interaction energies) may be unable to liberate them. Metastability and 
kinematic constraints ensure that dislocations are, and tend to remain, in the crystal despite the fact that in thermal 
equilibrium they should not be there. This raises major conceptual problems. Unless methods for dealing with 
metastability and constraints in statistical thermodynamics are developed to a much higher level than presently 
available, it may be difficult for us to derive the non-equilibrium statistical properties of dislocation systems from 
first principles. The present work therefore pursues more modest goals: To establish the fundamental structure of 
functionals which express the energy of a dislocation system in terms of dislocation density functions, to clearly 
formulate the parts of the functionals which we can know for certain, and to find reasonable approximations for 
those which depend on information regarding the relative arrangement of dislocations. In this task, we are inspired 
by density functional theory of electron systems where energy contributions of known and established form which 
represent long-range electrostatic interactions (the Hartree energy functional) are separated from those which depend 
on correlations in a many-body problem (the exchange-correlation energy), see e.g. 2^ or other textbooks on 
quantum mechanics of many-electron systems . The latter are approximated by reference to idealized systems such 
as a homogeneous electron gas. We apply the same strategy to the many-dislocation problem. Just as in electron 
theory, what we get is the fundamental mathematical structure of the energy functional which we determine first in 
two and then in three dimensions. As to the correlation energies, these depend on parameters characterizing the range 
and nature of dislocation-dislocation correlations - parameters which remain to be determined by reference to direct 
DDD simulations of the many-dislocation problem. Again, this strategy is analogous to the proceedings in density 
functional theory where exchange-correlation energy functionals are formulated and parameterized by reference to 
idealized model systems (thee free electron gas) or by comparison with direct numerical simulations of the many-body 
problem via quantum Monte Carlo methods. 

The technical method which we shall use is to represent the dislocation interaction energy in terms of densities and 
correlation functions. This idea is not new: Correlation functions have been introduced for averaging the forces in 
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2D dislocation systems, and hence the dynamics, in earlier work by Zaiser and Groma Averaging the forces 

implies, of course, evaluating average derivatives of the elastic energy functional. Here we apply the same averaging 
method to evaluate the energy functional itself. A similar approach has been used already in work the 1960s, see e.g 
the work of Kocks and Scattergood [ 2 ^ on systems of straight parallel dislocations, but was not further pursued. A 
possible reason for this lies in the fact that, in absence of DDD simulations which can provide complete information 
about the dislocation microstructure, information about dislocation correlations is hard to come by - even though 
some limited information can be inferred from electron microscopy and X-ray profile analysis data, and indeed some of 
the early work by Wilkens in the field considers both mean square stresses (and thus elastic energy densities) and X- 
ray line broadening (^ . . In the present study we resume these approaches and extend them to general dislocation 

systems in three dimensions with multiple slip systems and arbitrarily curved dislocations. We first revisit results of 
classical dislocation theory for the energy of discrete dislocation systems, and then develop our averaging methodology 
for the conceptually simple case of systems of straight parallel edge dislocations. We generalize the results first to 
curved dislocations on a single slip system, and then to general 3D dislocation systems. We then demonstrate a few 
applications for the resulting free energy functionals, first to evaluate ’Debye screening’ of dislocations in 2D, i.e., 
the formation of an induced distribution of local excess Burgers dislocations around a given dislocation which screens 
the long-range dislocation stress field. We also investigate the emergence of ’back stress’ terms in the stress balance 
that are proportional to second-order in the plastic strain gradients and demonstrate that these are associated with 
energy contributions that are quadratic functionals of the local excess dislocation density. Finally, we outline how our 
results can be used to evaluate the ’friction stress’ associated with formation and breaking of junctions in general 3D 
dislocation systems. We conclude with a brief discussion which puts our results into the context of other published 
work. 


ENERGY OF A DISCRETE DISLOCATION SYSTEM 


As demonstrated by de Wit [8| , the energy of a three dimensional system of dislocations can be written in terms of 
double integrals over the dislocation lines. This representation has been directly implemented in discrete dislocation 
dynamics (DDD) codes, notably the parametric DDD model of Ghoniem and co-workers 28, 2^ who use the variation 
of the energy with respect to dislocation line parameters in order to derive generalized forces acting on the dislocation 
lines. In our presentation we follow the representation given in the standard textbook of Hirth and Lothe Q which 
can be applied both to closed loops and to loop segments. We write the energy of a dislocation system consisting of 
closed loops with Burgers vectors as 


2 ^ /c(o /co) 


( 1 ) 


Here the sum runs over all pairs of loops, and the self-energy of each loop is evaluated as half the interaction energy 
of two loops at distance b (more generally, the core radius). Here and in the following, upper bracketed indices 
enumerate dislocation loops or segments of loops, whereas lower indices indicate coordinates of a Gartesian coordinate 
system. When dealing with 3D dislocation networks, we retain the decomposition into closed loops but break the 
loops into segments S separated by nodes. A junction which forms at the intersection of two loops and is 
thus represented as two segments of Burgers vectors 6^*^ and that are aligned with each other between the two 
nodes which delimit the junction. A collinear reaction where b*-*^ = —is represented as two aligned segments of 
opposite Burgers vector, not as a missing segment. In segment representation, the energy of the dislocation system is 

E=ly" f f -r(j)).Z(^')dsWds(^). (2) 


In Eqs. © and the interaction kernel is given by 
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In this expression, ^ is the shear modulus, v is Poisson’s number, G is a tensor with components 
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and has the components 
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In equation (|3]), the first term on the right-hand side is non-zero only if neither the line directions nor the Burgers 
vectors of both segments are aligned with each other, hence it can be considered to describe edge-screw interactions. 
The second term on the right-hand side describes the interactions of the screw components of both line segments, and 
the third term describes the interactions of edge components. 

These equations apply to dislocations in an infinite medium. In the presence of boundaries and boundary tractions 
which cause, in a fictitious crystal without dislocations, the stress field (Texti'f'), the energy changes. The corresponding 
energy contribution can be written in terms of the fictitious work that would need to be done by the Peach-Koehler 
forces in order to expand the loops C to their current size. For a system of planar glide loops this is simply given by 

Ee.t =y, [ : (Text(r)d2r. (7) 


Here is the area enclosed by the loop in the slip plane with normal and the projection tensor = 

0 + 'nA'> 0 e^'' ^)/2. Alternatively, we may write the same expression in terms of the microscopic plastic strain 


-pi,d 


(r) = 5{r — r')d^r'. 

4 Ja(‘) 


( 8 ) 


Inserting into Eq. © gives 


E'ext = [ eP'’‘^(r) : crext(r')d^r-. (9) 

Jv 

where the integration is carried over the crystal volume. In the following it will be useful to develop a number of ideas 
first for the physically unrealistic, but conceptually simple case of quasi-two-dimensional (2D) systems consisting of 
straight parallel edge dislocations pertaining to a single slip system. For such a dislocation system, we may without 
loss of generality set and where € [1) ~1] is the sign of a dislocation. The line integrals then 

reduce to integrals over the z axis, and the energy of the system becomes 

E = +y]F;,eif (10) 

i^j i 


where the vectors rO) now lie in the xy plane and all enemies are understood as energies per unit length in ^ direction. 
The self and interaction energies are given by (see e.g. [9|) 


Eself 
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( 11 ) 


where R is the crystal radius (or another external dimension of the system) and r = {x^ -f is the spacing of 

the dislocations in the xy plane, ab is a measure of the dislocation core radius, with the parameter a 1 chosen 
to correctly represent the core energy of the dislocations. Finally, the external energy for the considered dislocation 
system reads 

Text(a;,y^*^)da; (12) 



where Text = M : cText with the projection tensor —{bx 0 -b 0 ea;)/2 is the resolved shear stress in the slip 
system. In Eq. (HU the glide distance travelled by dislocation Q) ranges from an arbitrary reference position Xq to 
its current position x^i). 
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DENSITY FUNCTIONAL THEORY OF TWO-DIMENSIONAL DISLOCATION SYSTEMS 


To write the energy of our 2D model system as a functional of the dislocation densities, we define discrete densities 
of dislocations of sign s, and discrete pair densities of dislocation pairs of signs (s, s'), as 


P^ir)= 

—s 


/Os),?(r, r') = Y, S{r - . 

j-gU) —g 


(13) 


With these, the energy of the dislocation system can be written as 

E = YJ = Eseifd^r+^Y^^' JJ Eint{r - r')p'l’^^{r,r')d^rd^r' (14) 

s ss' 

where N is the total number of dislocations. We now make a transition towards continuous densities via an averaging 
operation (...) (see Appendix A). This leads to 

E = Y J Ps{i')Ese\{d‘^r+^Y^^'J J Eint{r - r')pss'{r,r')d‘^rd'^r' . (15) 

s ss' 

Here, the averaged single-particle densities Ps(r) = (Ps(r)) can be understood as averages of the sign-dependent 
discrete densities, and the averaged pair densities pss'{r,r') = {pi{r)p'^,{r')) are averages of products of discrete 
densities. Note that the averaged pair densities are in general not equal to the products of the averaged single¬ 
dislocation densities: averaging is a linear operation which does not interchange with the formation of a product. 

Hence, the information contained in the single-particle densities is of necessity incomplete. Nevertheless it is our 
goal to express the energy functional in terms of the densities Ps{i')- To this end we write the pair densities without 
loss of generality as 


Pss'{r,r') = ps{r)p,,{r')[l + dss'{r,r')] , (16) 

where dss> are correlation functions. This allows us to split the energy functional into a part which can be exactly 
expressed in terms of the dislocation densities (’Hartree Energy’ Eh), and a part which depends on the correlation 
functions (’Correlation Energy’ Eq) and needs to be evaluated in an approximate manner. After some algebra we 
arrive at 


E = Es + Eh + Eg 

= J p(r)Aseifd^r -t- ^ JJ K{r)K{r')Eint{r - r')d^rd^r' 

+ i^ss'^Ps(r)ps'(r’')4s'(’'!'^')^int(r-r')d^r-d^r'. (17) 

ss' 

Here we have introduced the notations 


Pir) = Y 

S S 


(18) 


for the total and excess dislocation densities. 


The Hartree or self-consistent energy 

The Hartree energy Eh depends only on the excess dislocation density. To analyze its meaning, we use that 
K = —{l/b)dx"f where 7 is the mesoscopically averaged shear strain on the single slip system. We may now integrate 
the expression for Eh twice by parts to write the Hartree energy as a functional of the mesoscopically averaged plastic 
strain eP = TWyP: 




eP(r) : r(r — r') : e'^{r')d^rd^r'. 


(19) 
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Here, for this particular problem, F = {M ^ ® M We can re-write Eq. (IT^ as 

Eh = y £^ir) : (Ti„t(r)d^r 


( 20 ) 


where the internal stress field is given by 


a-intir) = J T{r - r') : e'P{r')d^r' . ( 21 ) 

This corresponds to the solution of the elastic eigenstrain problem in an infinite medium by means of a Green’s 
function method, see e.g. [s^ where expressions for F are given for the case of a general plastic strain field. Thus, 
the Hartree energy is just the elastic energy associated with the elastic-plastic problem in the absence of boundary 
effects. In the general case where boundaries are present, the boundary conditions result in an additional, ’external’ 
stress field <Text which superimposes on the internal stress field cri„t and which enters into the energy Eext- In most 
practical circumstances where elastic-plastic problems are to be solved, the Hartree energy or its functional derivative 
(the internal stress) will not be computed from the dislocation field k but evaluated in conjunction with the ’external’ 
stress from the solution of the elastic boundary value problem. 

In summary, the Hartree Energy represents the part of the elastic energy functional that is related to long-range 
internal stresses which can be described in terms of the coarse grained strain field £P or its spatial derivative, the 
geometrically necessary dislocation (GND) density k. In absence of mesoscale strain gradients (k = 0, equal numbers 
of dislocations of both signs), this term is zero and, hence, all interaction energy terms are associated with correlations. 
We now focus on the correlation energy Eq ■ 


The correlation energy 

To evaluate the correlation energy we proceed in the spirit of density functional theory of electron systems, i.e., 
we use a local density approximation where we approximate the correlation energy by a local functional of the 
dislocation densities which we evaluate as the correlation energy of a spatially homogeneous reference system. As 
shown in Appendix B, this is feasible if and only if the correlation functions dgs')*’, f"') are short ranged functions of 
range i which, for large values of |r — r'| ^ £, go to zero faster than algebraically (short-range correlated/macro- 
disordered dislocation systems). The assumption that the dss' are short ranged functions allows us, for dislocation 
arrangements where the densities ps are weakly space dependent on scale i, to approximate 


Ps{ri)ps'{r2)dss'{ri,r2) Ps{ri)Ps'{ri)dss'{ri - r2). (22) 

This local density approximation represents the zeroth order of a systematic expansion which expresses the energy 
functional in terms of gradients of the dislocation densities of increasing order (see Appendix). From Eqs. (1171) and 
(1^ the correlation energy reads 


Ec = ^'^ss' J ps{r)ps'{r) J Eint{r)dss'{r)d'^r' 


d^r 


Using the notation 


this can be rewritten as 




Ec = 


27r(l — v) 


47r(l — z/) 




.ss'—zLl 


where 


Tss'= j ss'dsB'{r)g{r)d^r. 


(23) 


(24) 


(25) 


(26) 
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To proceed further we need to specify some properties of the functions dss' ■ By construction (see Appendix A) the 
correlation functions have the properties 


dss'ir - r')d^r' = 0, if s s' , dss{r - r')(rr' ^-1/ps{r) if s = s'. 


(27) 


where we have used that, in a weakly heterogeneous dislocation arrangement, Ps(r’) ~ Psi'f'') changes little over 
the range of the function dss- (This weak heterogeneity condition is a requirement for the use of a local density 
approximation, see the appendix for an outline towards more general non-local approaches). Furthermore, as discussed 
in detail elsewhere IJ, l3l| , due to the scale free nature of dislocation-dislocation interactions in nearly homogeneous 
dislocation systems, any correlation functions which emerge spontaneously from the evolution of an initially disordered 
dislocation system must exhibit a range i oc that is proportional to the mean dislocation spacing, i.e., to the 

inverse square root of the total dislocation density as defined by equation m- Hence we assume that the correlation 
functions depend only on the relative position of the two dislocations divided by the mean dislocation spacing, 
expressed through the variable u = (r — r')^, i.e. dss'{r — r') = dss'{u). Note that, since Ps{r) ~ Psir') over the 
range of the correlation function dss' , it does not matter whether we evaluate p at r or at r'. 

Let us now first consider the case of s ^ s'. As the first step the radial function g given by Eq. ((Ml) is rewritten as 


gir) = go + griryfp) where go = In {Ryfp) , = - In [vy/p) - . 


(28) 


After substituting Eq. (1281) into Eq. (1261) and using the condition (IMl) one obtains that 


Tss' = — 


= - d. 
P J 


dss'{r)gr{r)d'^r 

(u) 


?/ 

ln(M) -I- -| 


d^u = 


Ds 


(29) 


It can be seen that the correlation function enters our further considerations only in form of the dimensionless number 
Dss'. 

In the case s = s', we substitute Eq. (IMl) into Eq. (IMl) . Due to Eq. (EZD we find that 


Tss = J dss{r)gr{r)d^r - ^\n{pR^) 


I 


dss{u) 


— In(it) — 


d?u — -— In (pR^) 
^Ps 


p 2ps ^ 


(30) 


By substituting Eqs. (I29I30I) into Eq. (1251) we get 


Ec = 




4tt(1 — z/) 
47r(l — v) 


(D+. + D_+) 
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Dl 1 , / p2'i^ I D]\K^ 2 
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Dll 
4 p 


d?r, (31) 


where Di = X)ss' ss'Dss' and Du = J2ss' Dss'- Here we have used that p± = (p ± k)/2 and, in a weakly polarized 

dislocation arrangement, D+_|_ = D _ The non-dimensional parameter a = exp(ZIi/2) can be envisaged as a 

dislocation screening radius, measured in units of mean dislocation spacings. We see that the correlation energy 
in local density approximation is indeed a local functional of the dislocation density functions. It consists of two 
contributions: first, the term proportional to p can be envisaged as a screening energy which reduces the dislocation 
energy as compared to a random dislocation arrangement. Second, the term quadratic in k characterizes modifications 
to local screening as we move from an unpolarized to polarized dislocation arrangements. 

The total energy function of the dislocation arrangement then reads 


E = 


pb"^ 


47r(l — ly) 


-J / ^ 7 "^^ // '^ir)>iir')g{r-r')d^r' 


(32) 
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Here the normalization constant po^ p is given by po = j{a^b'^). The total energy consists of three contributions: 

The term proportional to p can be envisaged as line energy of a dislocation which is screened by the surrounding 
dislocations. The quadratic but local term in k is a correction to screening, and finally, the Hartree Energy expressed 
by the double integral over k describes the energy stored in the long-range elastic field associated with macroscopic 
polarization of the dislocation arrangement. 


DENSITY FUNCTIONAL THEORY OF THREE-DIMENSIONAL DISLOCATION SYSTEMS 

System of loops on a single slip system 

We first consider a 3D system of loops pertaining to a single slip system with slip vector e^^'> = and slip plane 
normal n = By. To facilitate the transfer of our results to the general case of multiple slip systems, we start from the 
‘segment representation’ of the dislocation energy, Eq. ©• For loops on a single slip system, there exists no natural 
subdivision into segments as provided by the nodes in a 3D dislocation network. Instead, the length of the segments S 
is an artificial parameter of the calculation, which will be chosen as a small fraction p of the local radius of curvature 
Rc of the dislocation lines - a method also used in discrete dislocation schemes which represent the dislocation as a 
sequence of straight segments separated by nodes [s^. We split the elastic energy of the system, Eq. © , into sums 
of segment self energies and segment interaction energies: 

E = Es + El 

1 


= / lis)-£ir{s)-r{s')).l{s')dsds' 

^ , Jsi') Js(') 


2 ^ 75 (*) Js(i) 


(33) 


Here, £ is the interaction energy tensor for segments in the considered slip system which we write as 

£{r-r') = !^g{r-r') 


(34) 


where g follows from Eq. (jT|) with ei = 62 = e^,. Since we consider glide loops on a single slip system where the line 
direction I is contained in the plane y = 0, this tensor has only two relevant components which are explicitly given by 


gxx{r - r') = 


1 


9zz{r-r') = — 


1 92 


V dy‘ 


r — r . 


(35) 


In a first approximation, the segment self energies are replaced by the self energies of straight segments which we 
write, following Hirth and Lothe Q, as 




(36) 


where the line energy tensor is given by 


F 


9l = 


(L ^B^x ® ^a;) 
1 - ^ 


In 


(37) 


In Eq. denotes a straight seg ment of length L connecting the endpoints of segment 5^*^. In the spirit of 

nodal discrete dislocation dynamics |32| , we take the segment length to be a small fraction of the radius of curvature 
of the dislocation line under consideration, L = rj/k where the dislocation curvature k characterizes the line shape 
on scales above the segment length. The deviation between the curved segment and its straight approximation goes 
to zero in proportion with 77 ^ as 77 —>■ 0. In practice, 77 may be adjusted to provide an optimum representation of the 
core energy contribution to the dislocation self energy. 

We now proceed in direct generalization of the two-dimensional case. We define the discrete dislocation segment 
density as 


Pd(T',(/') = V [ 6{r - r{s))6{cj) - cj){s))ds , 
c JS 


( 38 ) 
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where (j) is the angle between the Burgers vector b and the line direction 1. Similarly, we define the segment pair 
density as 

^ (( 5{r - r{s))5{r'- r{s'))5{<i) - - (j){s'))dsds' (39) 

Jss' 

This allows us to write the dislocation self energy as 

Es= JJ m.£LmPd{r,cl))d\dcj) (40) 

and the dislocation interaction energy as 

Ei = ^ JJJJ [l{(l)).£{r - r').l{(j)')] pP(r, r', (/»,())') dVd(/)dV'd(()' . (41) 

Upon averaging, the discrete densities and pair densities become continuous functions of their arguments. In Eq. 
(HOI) . the only change is that we replace the discrete density pd{r,(j)) by its continuous ensemble average p{r,(j)). To 
evaluate the dislocation interaction energy, we proceed in analogy with the 2D case and write the ensemble averaged 
pair density as 


Pp(r, (t), r', (j)') = p{r, 4>)p{r', (j)')[l + d{r, r', (f), (j)')] . (42) 

Comparison demonstrates that the angle 4> in the present formalism plays very much the same role as the ’sign ’s in the 
2D problem. Inserting Eq. (I42|) into Eq. dm allows us to separate the dislocation interaction energy Ei = Eh + Eq 
into Hartree and correlation energy terms. The self energy, Hartree energy and correlation energy are given by 

Es = JJ p{r,(j))[l{(l)).£L-l{(l))]d^rd(j) (43) 

Eh = ^ JJ p{r,(j))[l{(j)).£{r - r').l{(j)')] p{r',(/)') d(l)d(j}' d^rd^r' (44) 

Ec = ^JJ p{r,(l))[lW-^{'f'- r').l{cl)')] p{r',<j)')d{r,r',(j),(j)') dcl)d(j)' d^rd^r' (45) 


We can eliminate the explicit angular dependencies from these equations by choosing an appropriate representation 
of the dislocation density functions. To this end we resort to the concept of so-called alignment tensors. 


Alignment tensor representation of the disloeation density funetions 

To simplify the expressions it is convenient to use an idea of Hochrainer Q and represent the angular dependence 
of the dislocation density function p{r, (jf) in terms of an alignment tensor expansion. Eollowing Hochrainer, we define 
the sequence of (reducible) dislocation density alignment tensors as 


p[°l(r) = 

1 p{r,(l))d(j) =: p{r), 


II 

j p{r,(j))l{4>)d4) =: /«(r). 


p[2l(r.) = 

1 p{'r,4>)l{4>)®l{4>)d(l), 


II 

1 p{r,(l))l{(j))[^l{(j))r~^d(l). 

(46) 


where denotes a n-fold tensor product with the line direction vector 1. The zeroth-order alignment tensor 

is the conventional dislocation density (line length per unit volume). The first-order tensor (or dislocation density 
vector) has as its components the edge and screw contributions of the geometrically necessary dislocation density. 
The second-order tensor contains the information about the distribution of the total dislocation density over edge and 
screw orientations, and so on Q. 
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The dislocation density function p{r,<f) can be recovered from the alignment tensors as follows Q: We denote 
the irreducible part of the tensor as with components . Furthermore, we denote as pi"! {4>) the n-fold 

contraction of pi"! with the direction vector l{4>). Then, 


p{r, (r) + V 2"pl"l (r, , pN (r, 0) = p'"] {r)k, (</>)... (^) . 


(47) 


The double angular dependency of the pair correlation function d{r, r', cj), (j)') can be expressed in terms of a double 
alignment tensor expansion. We define 




(48) 


This is an expansion on the direct product of the unit circle with itself as the I are unit vectors in the xy plane. In 

— ?77,j 

terms of the associated irreducible tensors d ’ the function d can be represented as 

. OO 

d(r,r',</),</.') = ^ ^ 2"+'"JI"’™l(r,r',</.,</)') , (49) 

n,m—l 


with the expansion coefficients 




h') = 




(r, r%, (</))... li„ {(j)) Ij, {(j)')... Ij^ {(j)') . 


(50) 


The self energy 

The different contributions to the energy of a dislocation system can be expressed in a natural manner in terms of 
the dislocation density alignment tensors. As immediately seen from Eq. (H41l and the definition of the second order 
dislocation density alignment tensor, the self energy of a dislocation system can be represented as 

Es= J ■■ pl^!(r)d^r = ^ 

The term under the integral has the character of a local energy density which is evaluated as a double contraction 
of the line energy and second-order dislocation density alignment tensors. It depends both on the local density of 
dislocations and on their character (edge/screw). 


The Hartree energy 

According to Eq. (H41) and the definition of the first order dislocation density alignment tensor pl^l (also termed 
dislocation density vector k), the Hartree energy can be represented as 

Ah = - JJ pl^!(r).5(i—r').pl^!(r')d^rd^r' (52) 

Alternatively, we might express this energy in terms of the classical dislocation density tensor cx (the curl of the 
plastic distortion). To this end we note that ol = pl^! ® 6 = k® 6 and define the fourth-rank interaction energy tensor 
'K = b®E®b to write 


Ah ~ 2 11 ~ r').a'(r')d^rd^r' . 


This is the formulation used by Berdichevsky 


(53) 
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The correlation energy 


To evaluate the correlation energy Eq, we make the same crucial approximation as in the case of the 2D dislocation 
system: We use a local density approximation based upon the idea that the correlation function d is short ranged, and 
that we may approximate the correlation energy by the correlation energy of a homogeneous system. In lowest-order 
approximation we thus set 




p(r, (l)')d{r - r', (j), 4>') 

on+m 

p{r,cl>)p{r,(j)') (r-r', (/),(/)') . 

m,n—0 


(54) 


We now insert this approximation into the correlation energy, Eq. (1451) . The form of the expansion coefficients given 
by Eq. (1501) then leads us to define a sequence of coupling tensors with components 


rp[n,m] 


• • -jm 


47r^ 






(r)d^r. 


(55) 


where we again have assumed short-ranged correlation functions to ensure existence of the spatial integrals. We find 
that the correlation energy can be written in terms of these coupling tensors as 




Stt 


fib 

Stt 


n,m—l 
2 oo 


n,m—l ^ 


(m) 


(56) 


where : denotes a n-fold contraction. We note that, for reasons of symmetry, all interaction tensors must 

vanish where n -I- to is an odd number. Such tensors involve, in Eq. 1561 odd numbers of products of the director 1. 
Since I changes sign under coordinate inversion (r —>■ —r), whereas the energy does not, all terms in Eq. dlSD with 
odd numbers of products of I must be zero. 

In case where dislocation correlations emerge from the evolution of a otherwise scale free dislocation system we 
expect them to obey the relation d = d{r^) =: u. Eurthermore we note that g{r — r') = ^g{u — u'). Using these 
relations we can write the coupling tensors in analogy with Eq. (1291) as 


rji[n^rn\ 


jj[n,m\ 


P 




2’Ti+n—2 

47r^ 




(«)5in.Jm(“)d^'W- 


(57) 


where u = 


Energy functionals for continuum dislocation dynamics 

In practical terms, it is desirable to truncate the alignment tensor expansion at some low order. This is done in 
continuum dislocation dynamics theories which represent dislocation systems in terms of the evolution of dislocation 
densiW alignment tensors, restricting themselves to the alignment tensors of order zero and one @ or orders one and 
two [7|. By evaluating the interaction coefficients using data from discrete dislocation dynamics simulations, 

we might arrive at unbiased estimates to which degree such truncated expansions faithfully represent the energetics 
of dislocation systems. 

We explicitly give the energy functional for an expansion containing the first and second order dislocation alignment 

[2l [21 

tensors. The diagonal components of the second order alignment tensor ps := pxx and pe := Pyy correspond to the 
screw and edge dislocation densities. Furthermore we use that, because of invariance under the transformation 
{(j) —> the correlation alignment tensor is diagonal with only non-vanishing components 

and d^yy^K The same is true for the interaction energy tensor with the only non-vanishing components pxx and pyy. 






12 


With these notations we write 


E = Es + i?H + Eq 
^62 


Itt 


Psir) 


1 - 


In 


Pir) 

Poir) 


d^r + 


fib'^ 

Stt 


K{r).g{r — r').K{r')d^rd^r' 


+ 


f pIDss + plDee + DsePsPe ,3 , pb'^ f K{r).D^^’^\K{r) 


Stt 


p{r) 


d^r + 


Stt 


p{r) 


-d-^r. 


The interaction coefficients for the dislocation densities are given by 


£»s: 

De, 

Us, 


= Jga:a:{u)d'i^^^{u)tfu=^ j gxxiu)d{u,(j),(j)') cos (j) cos (j)'d?ud4>d(l)' 

= ^ J gzz{u)d^^^^\u)d^u = ^ J gzz{u)d{u,(j),(j)')sm4'sm(j)'d^ud(j)d4'' 

= X [[gzz{u)d'ii'^\u) + gxx{u)d^zi^\u)]d^u 


(58) 


d{u, (/>, (j>')[gyy{u) COS (j) COS (j)' + gxx{u) sin cj) sin (j)'](rud(j)d(j) 


(59) 


and the interaction matrix associated with the dislocation density vector = k is 




dTT^ 


g{u)d^°’°^{u)d^u = j g(tt)d(M, (/), (/)')d^ud(/)d^'. 


(60) 


The quantity po(^) = q{r)gb in Eq. (l5^ relates to the so-called curvature density (a product of dislocation density 
and curvature) which is one of the field variables of continuum dislocation dynamics as introduced in Q . 

A simplified theory as proposed in Q might only consider the dislocation density alignment tensors of order zero 
and one, hence retains only information about the total dislocation density p and geometrically necessary dislocation 
density vector k. Since no additional information is available, the alignment tensor of order two is then represented 
as = {p/2)I. With this simplification the energy functional becomes 


E = 




p{r )In 


pb^{2 — v) 
87r(l — v) 


Pb^ f /^(r).£)^"’"bK(r;_„3 
Stt j p{r) 


Pjr) 

Pair) 

-d“r 


pb"^ 


d'^r-I-— II K{r).g{r — r').K{r')d^rd^r' 


(61) 


Here, all coefficients of energy contributions proportional to p have been absorbed into the scaling factor pg = aqb, 
where the numerical parameter a is proportional to exp[(Ilss + -Dee + Dse)/4]. 


Multiple slip systems 


In case of multiple slip systems, dislocations are likely to form 3D networks. We now consider as segments 5*’^ 
stretches of dislocations bounded by two nodes in the dislocation network. The superscript /3 distinguishes the different 
slip systems with Burgers vectors be^, slip plane normal vectors and projection tensors . The interactions 
between segments pertaining to two slip systems /3 and /?' are given by where = e^, . For 

self energies we again approximate the self energy of a segment by that of a straight segment, for which we introduce 
the line energy tensor 



9l = 


l-iy 


{I - 0 e'^) In 



1 

2 ( 1 -^^) 


{I 


ve^ 0 e^) In 



(62) 


where we used that the segment length (mesh length of the dislocation network) is now proportional to the character¬ 
istic dislocation spacing (1/^/p). The parameter g can again be adjusted to account for the dislocation core energy. 
There exists, in a three dimensional network, the possibility that a segment of a loop of slip system /3 is collinear with 
a segment of slip system /3' (the two segments form a junction of Burgers vector ). We evaluate the 
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junction energy as the sum of the energies of the constituent segments and an interaction energy. This interaction 
energy is strictly negative (otherwise the junction does not form). It is given by 


E 


/3/3' 

J 


_ „/3/3' 




(63) 


where is the direction of the junction segment which is constrained to form along the line of intersection of the 
slip planes of both slip systems. At first glance, our method of book-keeping may look unusual - why not directly 
evaluate the self energy of the junction segment and discarding the addition-cum-subtraction of the constituent 
segment energies? The reasons for this procedure will become transparent later. 

Discrete dislocation densities are now defined separately for each slip system as 


Pdi'^A^) = - r{s))6{(l)^ - (l)^{s))ds 


(64) 


where is the angle between the Burgers vector and the line direction l{s). For junction segments we define 
junction densities 


p^^{r,(j)^') = 

Here, the function has the value 1 whenever a segment of slip system /? forms a junction with a segment of slip 
system /3', and the value 0 otherwise. Note that a junction can alternatively be envisaged as a segment of orientation 
(pPP = arccos(I^^ e^) in slip system /3 or as a segment of orientation cfi^ ^ = arccos(i^^ ) in slip system j3'. 

For describing interactions between non-collinear dislocation segments, we define pair densities for pairs of slip 
systems as 

{r, (j)^, r \^ ff 5{'i — r{s))5{r' — r{s'))5{4i^ — (j)^{s))5{(j)^ — (s'))dsds'. (66) 

Upon averaging, all these densities become continuous functions of their arguments and we drop the subscript d. As 
previously we write the pair density functions in terms of products of single dislocation densities and pair correlation 
functions, 

p^l^' = pf^{r,(j)^)p^{r,(j)f^)[l + ,r',(j)^')] (67) 

We skip the intermediate steps which proceed in direct analogy with those for a single slip system, with the only 
differences that now we need to sum over all slip systems (for the self energy) and all pairs of slip systems (for the 
interaction energy), and that we need to account explicitly for the junction energy. As in the previous section, we 
expand the slip system specific dislocation densities and correlation functions into alignment tensors 


p0A^](^r) = j p^{r,(j)^)l{4'^)[iSili4>^)]'^ ( 68 ) 

Using these notations we write the self and Hartree energies as 

(70) 

Ah = ^ [[ p^’^^^{r).g^^ {r — r').p^’^^\r')d^rd^r (71) 

/3/3' 


/ 6{r - r{s))S{(l)^ - (l)^'^')ds = p'^{r,cl)^)f>^>^'{r,cl)^) 

Y j 5{r - r{s))5{(j)^' - Y^{s))ds = p^^{r,Y)f^'^{r,Y) (65) 
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where is obtained from Eq. o by setting = e^, . A new contribution to the system energy in case 

of multiple slip systems is the junction energy. To represent the junction energy we introduce the definition 




fi' 


(72) 


with ^ <\ to meet the conditions 'Yhp' < 1 (only a fraction < 1 of all dislocation segments can form 

junctions) and ^p^ (a junction between segments of slip systems f3 and /3' is a junction of slip systems 

P' and /3). Using this notation we can write the junction energy as 


Ei 



pP{r)pP (r) 
p{r) 


/i^^'(r) 5 f'dV 


(73) 


In this expression we have made the simplifying assumption that the probability of forming a junction does not depend 
strongly on the orientation of the intersecting dislocations in their respective slip planes. A more general treatment 
which uses an alignment tensor expansion of , and of which Eq. (1731) is the lowest-order term, will be given 
elsewhere. 

We are left with evaluating the correlation energy which contains all terms dependent on the correlation functions 
. This can be written as 


PP' n,m=l'' 


where the interaction coefficients are 

rtPP'lrt^rn] 

f'>o •yon 1 


[n,mj 


2771+71 —2 

47r^ 


T/3/3'[n-l,m-l] 

'^1 ■ • •'^n — 1 :jl • • -jn 




(74) 


(75) 


with u = r^fp. We note that the main qualitative difference between the single and multiple slip situations resides in 
the possible existence of collinear segments, i.e. junctions. In comparison with mutual interactions between distant 
segments of different loops, junctions may lead to a much more efficient energy reduction. 


APPLICATIONS 


Dislocation screening in two dimensions 


As an application of our two-dimensional theory, we revisit the problem of Debye screening of dislocations which 
has been previously studied by Groma and co-workers [l^ . We use the energy functional given by Eq. (IH^ to evaluate 
the reponse of a homogeneous, infinitely extended 2D dislocation system of density p to a single excess dislocation 
fixed in the origin, ko(t’) = 5{r). The induced excess dislocation density k follows by considering the variation of the 
ensuing energy functional, under the assumption that the overall density p remains homogeneous. The variation of 
Eq. (15^ with respect to n is then given by 


SE = 


47r(l — I/) 


Dii^-k J[K{r') + S{r')]g{r - r')d'^r' 


SK{r)d'^r = 0 . 


This leads to the following equilibrium equation for the induced density k: 


pb'^ 

47 r(l — v) 


2 p 


K{r')g{r 


r')d^r' -I- g{r) 


= 0 . 


(76) 


(77) 


The solution of Eq. dZZl) can be found by Fourier transformation. Using 


ff(fc) 


Sirky 


(78) 
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and the definition = Anp/Du, we find 


K{k) = — 


Aklkl 


+ k^ 


from which reverse Fourier transformation yields the result 

2 /sinh(fcoj/) 


c(r) = 


^0 


Ki{kor) - cosh{koy)Ko{kor) 


(79) 


(80) 


This is also the result obtained by Groma and co-workers [lfil| and, using a quite different formalism, by Linkumnerd 
and Van der Giessen [^. We point out that our investigation, though it leads to the same result, differs somewhat 
from the work of Groma and also of Limkunmerd and Van der Giessen. Groma et. al use a free energy functional 
which is devised heuristically and the term D, which controls the range of correlations, is associated with entropy-like 
terms in the free energy. The same is true for the investigation of Limkumnerd and Van der Giessen 21| who relate 
the range of correlations to fluctuation terms in the dislocation dynamics which they characterize by an effective 
temperature. In the present investigation, on the other hand, the parameter /cq, or Du, which controls the interaction 
range, arises from purely energetic considerations. Comparison with discrete simulations gives fco = 4.2y^ [l^ which 
allows us to obtain the numerical value of the coupling parameter Du Ri 0.84. Inserting this numerical value into Eq. 
(EH) together with typical values p lO^^m ^ and po ~ h ^ ~ 10 ^®m 2 indicates that, in the absence of long-range 
stresses (Hartree energy), the composition of the dislocation arrangement (geometrically necessary vs. statistically 
stored dislocations) has only a quite modest influence on the energetics. Even in the extreme limit n = p (only 
geometrically necessary dislocations) the additional energy cost implicit in the term proportional to amounts only 
to about 3% of the term proportional to p. 


Derivation of ’back stress’ terms in two and three dimensions 


In dislocation-based plasticity theories, many authors have found it convenient to introduce ’back stress’ terms 

iilllF 


proportional to the gradient of the dislocation density vector k into the stress balance, see e.g. 


Such 


terms are of interest also because the dislocation density vector k is proportional to the gradient of plastic strain, 
hence, ’back stress’ terms correspond to second-order plastic strain gradients entering the stress balance, a device 
highly popular in phenomenological gradient plasticity models of continuum mechanics. We demonstrate in this 
section that such terms arise naturally from our density functional representation of the dislocation energy. 

To this end, we first consider the 2D case. We take the /t-dependent terms in the energy functional given by Eq. 
(EH) and insert the relation between the excess dislocation density k and the plastic strain 7 , k = —{l/b)dx'Y- 


E = 


167 r(I — i/)p 

Variation with respect to 7 yields 


/ ^ 8Tr{l- v) J - r')d'^rd^r' . 


(81) 


SE = 


47 r(l — v)p 


- J dxi{r')dxg{r - r')dV' 


dyd^r = / T{r)5^{r)d?r 


(82) 


where we have used that the work conjugate of the plastic shear strain 7 is a resolved shear stress in the considered 
slip system. Using that the shear stress of a single dislocation is given by rd(r) = —^ 6 /( 47 r(l — v))dxg{r) we see that 
the shear stress in the slip system is of the form 

t(’') = + f - r')d^r' = + t^c (83) 

87 r(l - v)p J 

The first of these terms is the back stress Tb derived, along a quite different line of reasoning, by Groma et. al [1^. The 
present derivation makes it obvious that this term results from the correlation energy contribution that is quadratic 
in the excess dislocation density k. The second term, t^c, represents the superposition of the long-range stress fields 
of the excess dislocations. This stress contribution derives from the Hartree energy and is normally obtained from 
solving the standard elastic-plastic problem. 

We can repeat the same argument for 3D systems. In case of a single slip system, the dislocation density vector 
relates to the strain gradient by k = (l/ 6 )e„V 7 , = e.n where e is the Levi-Civita tensor and n the slip plane 
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normal. The tensor e„ rotates a vector contained in the slip plane, such as k, counter-clockwise by 90°. With this 
notation we can write the K-dependent terms in the energy functional, Eq. (I58L as 


E(k) = - / / K(r).g(r — r ).K(r )d°rd°r -|-/ ---^d°r 

Sir JJ 8 tt J p[r) 




p{r) 


= ^ JJ[er,.V'r{r')].g{r - r').[er,.V'r{r)]d^rd^r' + ^ J 

Variation with respect to S'y gives 


p{r) 


(57(r)d^ 


The term in the brackets can again be understood as the resolved shear stress. In terms of n it is given by 


T(r) = [ K{r')[enV]rgir - r')d^r'- ^[enV].D^^’^\K{r) = Ti,+' 

Att j Airp 


(84) 


(85) 


( 86 ) 


In case of multiple slip systems, we use the notation = e.n^ where is the slip plane normal of slip system /3. 
After repeating the steps as above we get for the slip system specihc back stress terms in multiple slip conditions 

= -^[e£v].[^r>^^'[i'iU^'(r)]. (87) 

p, 


Some implications of our derivation of the back stress term are discussed in Appendix C. 


Estimate of the friction stress for a dislocation moving in a multiple slip environment 

Because of the geometrical constraints to dislocation glide on slip planes, dislocations can in general not move 
without intersecting dislocations on other slip systems. Because of this, sustained dislocation motion requires the 
repeated formation and breaking of junctions. The work required to break junctions is dissipated in the process. The 
distance between statistically equivalent configurations in the average direction of dislocation motion is given by the 
mesh length Ij^/p of the dislocation network. Thus, the energy dissipated in advancing the dislocation by a distance 
5u > 1/y/p can be estimated as 

^diss('5«) = X! X! f [Vpi5M(s)]ds =-J2f (.r)bSu{s)ds (88) 

where is the friction stress required to move the dislocations, i.e., the resolved shear stress r which provides the 
work required for breaking junctions. Upon transition to an averaged formulation we can re-write this equation as 

^Ls = 47,-(1 - ,y) ]Su{r)d^r = j p^{r)T^{r)6u{r)d^r. (89) 

Since the virtual displacement Su is arbitrary, it follows that the shear stress (’friction stress’) required to move the 
dislocation is 

"■f" y) ^ (90) 

This stress obeys the generic Taylor scaling relation, i.e., it is proportional to the square root of dislocation density. 
It also depends on the distribution of dislocations over the various slip systems and on the coefficients which 

have the character of latent hardening coefficients. The basic idea underlying the above argument is that the junction 
energy defines the amplitude of the small-scale energy fluctuations (on scales of the order of one dislocation spacing) 
that need to be overcome in order to move a dislocation by repeated breaking and formation of junctions. We note 
that the above argument can be generalized by replacing the scalars and p^ with alignment tensor expansions of 
the corresponding angle-dependent functions. In this manner one can account for the fact that the average junction 
length may depend on the orientation distribution of the intersecting dislocations. This will be discussed in detail 
elsewhere. 
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DISCUSSION AND CONCLUSIONS 


It is interesting to compare our results with related work by other researchers, notably regarding the structure of 
the energy functional. We have shown that the energy functionals of dislocation systems possess a generic structure 
which is common to 2D and 3D dislocation systems. Specifically, the energy functionals consist of a ’Hartree’ energy 
which is a non-local, quadratic functional of the dislocation density vector, or equivalently of the dislocation density 
tensor. This part of the energy functional does not depend on assumptions regarding dislocation correlations. The 
Hartree energy is complemented by an energy term which has the form Eg (x — J p\n{p/pQ)d^r where po oc (1/6^). 
This energy is proportional to the line length per unit volume with a proportionality factor that decreases with 
increasing dislocation density, reflecting the fact that the screening radius of dislocation systems is proportional to 
the dislocation spacing. Terms of the form —plnp in a free energy density are normally associated with entropy, 
and indeed such terms appear in thermodynamic theories of dislocation systems, see e.g. Kooiman [l^ . However, in 
thermodynamic theories the pre-factor of pinp type entropy terms is bound to be of the order of kT, which is several 
orders of magnitude less than the actual pre-factor Ri p,b^. The present derivation makes it clear that this term is in 
fact of energetic origin. We note that, in the hypothetical case where the dislocation density approaches po, according 
to the present formalism the energy per unit dislocation length in a system that is equally composed of positive and 
negative dislocations would go to zero. This is simply a reflection of the fact that in this case the cores of the positive 
and negative dislocations overlap and the dislocations annihilate. The role of the parameter Po is thus the exact 
opposite of the ’limit dislocation density’ ps introduced in an ad-hoc manner by Berdichevsky |37| . This term was 
introduced into the logarithmic factor in such a manner that it makes the energy per unit dislocation length diverge 
as the dislocation density approaches the critical value pa- In view of our results this idea must be discarded. Indeed, 
if we consider dislocation systems of zero net Burgers vector, it is difficult to see how densification of the dislocation 
system, i.e. bringing dislocations of positive and negative sign closer to each other, could conceivably increase rather 
than decrease the energy per dislocation length. The third energy contribution which consistently emerges from the 
present treatment is a local term which is quadratic in the excess (geometrically necessary) dislocation density. This 
term forms part of the ’correlation energy’; it depends on the structure of the dislocation pair correlation functions. 
Upon variation, it yields the ’back stress’ which has become very popular in both phenomenological and dislocation 
density-based plasticity theories, not least because of its ability to explain size effects [Il,|3^. To summarize, we hold 
the following fundamental structure of the energy density in the dislocation energy functional to be generic: 

• Self energy terms of form oc —J pln(p/po)d^r 


• A non-local Hartree energy which depends on the excess (GND) dislocation density vector, or equivalently on the 
dislocation density tensor, of form pb^ ff K(r).g(r — r')K(r')d^rd^r' or pb^ ff a(r).7Z(r — r')a(r')d^rd^r' 

• Terms proportional to the square of the excess dislocation density vector, of form pb^ f l^K.D.Kd^r 

Of these terms, the self energy depends only logarithmically (through the term po) on the correlation functions. Also 
dependent on the structure of the correlation functions are the length scale i and the coupling tensors D. Of these 
we only know that they must be positively definite, since otherwise a homogeneous system of statistically stored 
dislocations would spontaneously decompose - which it does not. As to the length scale i, in the absence of other 
factors it must, in order to be consistent with the scaling properties of discrete dislocation systems 13, 3lj, be chosen 
proportional to the dislocation spacing. This is the approach used in the present work. 

In our evaluation of the energy functional of a dislocation system we have made the key assumption that the range 
of correlations between dislocations is limited. This is a necessary assumption for expressing the correlation energy 
as a local functional of the dislocation densities, see Appendix B. In physical terms this assumption corresponds to 
the simple idea that, by studying the dislocation configuration in one point, one cannot gain any information about 
the configuration of dislocations in a distant point (at a distance of many dislocation spacings) which is not contained 
in the slip-system dependent dislocation densities in that point. While this seems a rather mild assumption, it must 
be noted that a lot of dislocation systems that have been extensively investigated in the literature fall not into this 
category, among them: 

• the Taylor lattice [s^. 


• the infinite periodic or non-periodic dislocation wall [35 

• the infinite dislocation pile-up 4fll. 411. 


• the periodic misfit dislocation array. 
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Many of these systems are one-dimensional and/or periodic, which makes powerful mathematical tools available for 
their analysis. The existence of these powerful tools may motivate the investigation but, on the other hand, it is also 
clear that most of the dislocation arrangements which develop during plastic deformation and whose properties govern 
plastic flow are disordered rather than ordered on large scales, and are two- or three-dimensional rather than one¬ 
dimensional. No Taylor lattice and no infinite periodic dislocation wall has ever been seen in the electron microscope, 
extended pile ups are the exception rather than the rule in deformation of real materials, and even in case of interface 
dislocations the assumption of periodic order has recently been called into question. We may thus argue that the study 
of low dimensional and/or periodic dislocation arrangements (which can never be captured by the present approach) 
is a consequence of mathematical convenience rather than of their practical importance. In this sense the present 
investigation can hopefully be considered a step in ’becoming generic’ which also means ’becoming realistic’. 

Among the applications we have given, we consider the derivation of back stress terms to be of fundamental interest. 
It should be clear from our derivation that the fundamental term in the energy functional that depends on the excess 
dislocation density k is the Hartree energy. Whatever assumptions are made regarding the correlation functions, 
this term is bound to stay. The term proportional to k? in the correlation energy, which gives rise to back stress 
terms, is in fact a local correction to the fundamentally non-local functional - a fact well recognized in recent work by 
Kooimans et. al. Accordingly, the back stress is a local correction to the in general non-local, long ranged interaction 
between excess dislocation densities in different parts of the crystal. In view of this fact it is astonishing that there 
are several published attempts to replace, rather than correct, the long-range dislocation interaction by a back stress 
term. This is tantamount to throwing out the Hartree energy and expressing the elastic energy as a functional of the 
excess dislocation density which does not contain any non-local, long-range interaction terms. The standard device 
used to this end is a truncation of the kernel g at some arbitrary radius £, see e.g. 34) We cannot help pointing 


out that this is inconsistent with the most fundamental property of dislocations, and of dislocation systems, namely 
the existence of a Burgers vector that is independent on the Burgers circuit: Let us compute the stress (or equivalently 
the elastic strain) associated with an arbitrary dislocation arrangement contained in some circle Cr of radius R. If we 
now truncate the stress around dislocations at length we are left with a horrible dilemma - either we truncate the 
stress but not the strain, in which case we have destroyed elasticity, or we truncate both, in which case the integral 
over the circle Cr+i will yield zero whatever the Burgers vector content in C might be, which would be possible only 
if dislocations had no Burgers vector to begin with. This idea is taken to its logical conclusion in the work of Luscher 
et. al. [ 3 ^ who use the back stress to evaluate the dislocation-associated strain as a compatible tensor field, see 
Appendix C. 

There are several directions how the present investigation could be expanded and further developed. At present, 
our treatment of segment self-interactions via a line energy approximation is not very elegant. This can be easily 
improved upon by replacing the interaction tensors g by core-regularized expressions which can be derived in various 
manners including gradient elasticty [l^ and continuous Burgers vector distributions around the dislocation core 
m- Since we cannot, in the general case, calculate dislocation correlation functions from our theory, we need to 
obtain the information regarding quantities like po and from external sources. In our opinion, as a next step, a 

systematic effort is needed to evaluate the expansion parameters of the theory (the coupling tensors D [n,m]-^ from DDD 
simulations. From a numerical point of view this is not difficult, especially with reference to DDD codes that express 
the interaction energy in terms of line integrals over the dislocation lines [ 2 ^, [ 2 ^ , since the coupling tensors can for 
discrete dislocation systems be evaluated in a similar manner. A comparison with DDD simulation data will serve two 
important objectives. Firstly, owing to the non-ergodicity of dislocation dynamics, it is not at all clear to which extent 
these coupling tensors depend on initial conditions. If they do so in a sensitive manner, meaning that different types 
of initial conditions lead to quite different values for the coupling tensors and thus to different energy functionals, 
then the present theory is useless for practical application. If, on the other hand, the dependence on initial conditions 
and deformation geometry is weak and remains within the statistical scatter among individual simulations, then the 
theory can be applied for evaluating the dynamics of dislocation systems from density-based evolution equations 
in a correspondingly wide range of situations. Secondly, comparison with DDD simulation can tell us how many 
terms of the alignment tensor expansion are actually needed for a meaningful representation of dislocation energetics, 
and thus provide important hints regarding the question what degree of complexity is actually needed for density 
based dislocation dynamics models. DDD simulations can be usefully complemented by experimental data regarding 
the structure of energy functionals for dislocation systems. Classical X-ray and calorimetry studies (see e.g. [H) 
provide information about the energy stored in a dislocated crystal and, via line profile analysis, about the rang e 
and dislocation density dependence of screening correlations in dislocation systems (Wilkens’ M-Parameter, [23j). 
Recent developments in X-ray microscopy allow to map the lattice distortions, and hence the elastic energy density 
associated with dislocation systems on scales well below the spacing of individual dislocations, see the impressive work 
of Wilkinson and co-workers ^ . Such experiments allow to obtain data that are of comparable quality to those from 
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DDD simulation and can be used in a similar manner for evaluating the elastic energy functional and parameterizing 
its local and non-local terms. 
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Appendix 


Averaging procedures 


To represent a discrete dislocation system by continuous densities, averaging procedures are required. Spatial 
averaging is commonly used in mechanics, but has the disadvantage that it does not preserve information about the 
relative positions of dislocations (or particles) with respect to each other, information which is essential for capturing 
the energetics. Hence, spatial averages are not normally used in statistical mechanics where one aims at deriving 
the average properties of systems from the dynamics and interactions of their discrete elements. The only averages 
normally used in statistical mechanics are (i) temporal averages along the trajectory of a system or (ii) instantaneous 
averages over ensembles of many systems. In thermal equilibrium, both types of averages are assumed to coincide 
(ergodicity). However, dislocation systems during plastic deformation are not in thermal equilibrium, and moreover 
their dynamics is strongly influenced by constraints (glide on crystallographic planes) which normally prevent them 
from fully exploring the phase space of possible configurations. Hence, dislocation motion tends to be non-ergodic. 
This leaves us with ensemble averaging as the only feasible averaging approach. 

Initial conditions for a 2D dislocation dynamics simulation are provided by assigning initial positions and signs to 
N dislocations. Any statistical rule for doing so explicitly or implicitly defines an initial N-particle density probability 
density function 

PAr(ri ... Tat, Si ... Sat, Si ... SAr)d'°ri . .. d-^rAr (91) 


which is the joint probability to find the first dislocation of sign 1 at ri, dislocation 2 of sign 2 at r 2 , etc. Obviously, 
the N-particle density function fulfils the normalization condition 

f p{ri.. .rN,si.. .SN)d^ri.. .d^rN = 1 . (92) 

SI...SN 


An non-equilibrium ensemble is defined by its initial probability density function and evolution equations, hence, 
the rules for constructing initial conditions in a set of multiple DDD simulations can be considered to define an 
ensemble. This is true for both 2D and 3D simulations: Initial conditions for a 3D dislocation dynamics simulation 
can be understood as statistical rules for assigning initial positions and directions to N segments in terms of a density 
function p{ri .. where those rules need to respect line connectivity and one makes the transition 

N ^ oo as the segments are made to be arbitrarily short. 

Probabilities of lower order can be obtained by integrating over some of the coordinates. Of particular importance 
in the present study are the single-particle and pair probabilities defined by 

pi(r,s) = pN{r,r 2 ■ ■ ■r^r,s,S 2 ■ ■ ■SN)d^r 2 ■. 

P2{r,sy,s') = J2 /Mr,r',r 3 ...r^,s,s',S 3 ...s^)d^r 3 ...d^r^. (93) 

S3...SAT 


These give, respectively, the joint probability of finding a dislocation at r and with sign s, and the joint probability 
of finding a dislocation pair at (r, r') with signs (s, s'), irrespective of the positions and signs of all other dislocations. 
The probability for a dislocation at any position to have sign s is given by p{s) = Ng/N where Ng is the number of 
dislocations of sign s. We write 

pi(r,s) =p(s)/f(r) , p 2 {r,r\s,s') = p{s)p{s')fl‘‘'{r,r') (94) 

where f({r) is the conditional probability density for a dislocation of sign s to be at r and /|® (r, r') is the conditional 
probability density for a dislocation pair of signs (s,s') to be at the positions {r,r'). From the general normaliza¬ 
tion condition, Eq. (1^ . we see that these conditional probabilities are normalized according to f ff(r)d^r = 1, 
fffi"'(r,r')d^rd^r' = 1 . 

From the sign-conditional single-dislocation and pair probability densities we obtain the respective dislocation 
densities by 


Ps(r) = iV^/f^(r) 


Ps<s'(r,r') 


NgNg, fy(r,r') , s'^ s 
iV,(iV,-l)/r(r,r') , s' 


(95) 
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Note that, while the number of pairs of dislocations of types s ^ s' is NgNs', the number of pairs of dislocations of 
type s is Ns{Ns — 1) since a dislocation cannot form a pair with itself. The densities are thus normalized to yield, 
upon spatial integration, the respective numbers of dislocations or dislocation pairs. From this normalization it follows 
that, if we express the pair density as Pss'{r, r') = ps(r)ps/(r')[l + dss'{r, r')] then 


J ps{r)dss'{r,r')d^r = J ps'{r')d,s'{r,r')d?r' 


0 , s' 

-1 , s = s'. 


(96) 


We see that the quantity ps{r)dssir, r') has a role similar to the exchange-correlation hole density in density functional 
theories of electron systems (see e.g. 2^). In a local density approximation where Ps{r) depends only weakly on r 
over the range of the correlation function dssir,r'), we may pull the factors Ps(r) ~ Psi'i'') out of the integrals, and 
Eq. (l27ll follows. 

In conclusion, a point about the choice of initial conditions (the intial many-dislocation density function, or the initial 
conditions in a series of DDD simulations) is appropriate. To enable comparison with experiment, initial conditions 
in DDD should be consistent with information about dislocation microstructure that is accessible by experiment. In 
well characterized microstructures, this is typically the total dislocation density, the geometrically necessary density 
as monitored by lattice rotations or misorientations, and possibly the distribution of dislocations over the various 
Burgers vectors. In extremely well characterized specimens, even information about the distribution of dislocations 
over edge and screw orienations may be available. 

All these informations are comprised in the dislocation density alignment tensors up to order two. However, any 
DDD simulation involves constructing initial conditions which imply the definition of a many-dislocation or many- 
segment density function - a function which may contain much more information than is contained in the simple density 
functions. It is in the opinion of the present author essential that initial conditions are constructed in a manner that 
does not introduce more information than is actually available, since otherwise the results may be influenced in a 
significant, and potentially uncontrollable, manner by hidden parameters introduced in the form of assumptions about 
the initial state that are not backed up by experimental evidence. A systematic manner of constructing unbiased 
initial conditions is provided by the maximum entropy method, which allows to construct the many-particle density 
by maximizing the entropy while using the available information as constraints. As an example, for a 2D dislocation 
system of size with A^-l- positive and N— negative dislocations (densities p± -\- /Li^), the N-particle probability 

density function which maximizes the entropy is simply pM = (I/T^)'^^(I/T^)^“ - in simple words, as an initial 
condition, the dislocations are placed independently at random locations, which is indeed a popular initial state 
for 2D DDD simulations. Other micro-arrangements which are popular in the literature for analysing properties of 
dislocation systems, for instance placing dislocations on regularly spaced slip planes [s^, or even on a regular Taylor 
lattic [ 3 ^, seem highly problematic from an information-theoretical of view because they imply strong assumptions 
about a correlation structure which may not be backed up by experimental evidence. 


Non-local density fnnctional approximations of the correlation energy 


The local density approximation used in this work can be considered the lowest order of a systematic expansion of 
the energy functional in terms of gradients of the dislocation densities. We illustrate this for the correlation energy 
of a 2D dislocation system. We start from Eq. (HZl): 

= J J Ps{r)ps'{r')dss'{r,r')Ei^t{r - r')d?rd?r' . (97) 

ss' 

We introduce the vectors r* = (r -t- r')/2 and a = 1 — r' and expand both p and ps' around the point r*: 


00 - 


n=0 


a.Vr 


Ps{r)\i 




m—0 




Inserting into the correlation energy gives 

2"+-ri!m! JJ (^ss'{r*,a)Eint{a) 

a a' r) m 


a.Vr 


Ps{r*) 


a.Vr 


Ps'{r)\r-. 


Ps'{r*)d^ad^r* . 


(98) 


(99) 
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We then introduce the gradient coefficient tensors T 


(n+m) 

ss' 


with components 


rpn+m /^*\ _ 

ss'.n.-.ir. + m'' > 2”+™n!TO! 


• ■ • Uj 


^dss'(r*, a)Eintia)d^a 


to write the correlation energy as 


( 100 ) 


^c = ^EE/ VV(r)^Tl:+™)(r)'^V>.,(r)dV. (101) 

ss' n,m 

Here, the tensorial m-th order dislocation density gradient V^ps is the rank-m-Tensor with components di-^ ... di^ps- 
Thus the correlation energy can be represented in terms of a gradient expansion of the dislocation densities, provided 
that the dislocation density functions can be differentiated to arbitrary order and that the gradient coefficient tensors 
of arbitrary order exist. A necessary and sufficient condition for this is a faster than algebraic decay of the correlation 
functions dgs'l^yV')^ which corresponds to the assumption of a macro-disordered dislocation arrangement. 

It is straightforward to generalize the above argument to 3D dislocation systems, however, the notation associated 
with a double expansion in real space and in angular coordinates is cumbersome so we refrain from giving explicit 
expressions. The local density approximation used in the remainder of this paper is just the lowest-order term of the 
above mentioned gradient expansion. The above argument demonstrates that this approximation can be systematically 
generalized to derive gradient-dependent expressions for the correlation energy of any desired order. 


How not to understand back stresses 


Eq. [ST] relates the back stress on a slip system to directional derivatives of the dislocation density vector: 


( 102 ) 


Our derivation tells us that this stress enters into the stress balance alongside the standard non-local stress (’Hartree’ 
stress). However, several authors [3J-|36| have suggested to use the back stress in order to replace long-range dislocation 
interactions. To illustrate the implications, we follow Luscher et. al. who use an expression exactly analogous 
to the above equation, though with a scalar coupling constant D. If this expression is assumed to fully describe 
the dislocation associated stress field, it is only natural to associate the back stress with a matching stress tensor, 
T? = and this with a strain through 


Cb = C-i : ^Tb = C-i : (M^)-Vb (103) 

where C is Hooke’s tensor. Luscher et. al. use this expression to evaluate the dislocation associated strain which 
is only logical, since there is no other stress associated with dislocations in their theory. To show the implications 
we look at a special case - our 2D dislocation system with Burgers vector b = be^ and slip plane normal n = By 
containing straight parallel dislocations of line direction I = e^. The dislocation density vector is k = kBz, and 
the back stress Tb = —{pb'^D)/{4Trp)dxK. We now consider a particular dislocation distribution - a blob of positive 
dislocations with density K{r) = kq exp(—r^/£^) and ask what is the Burgers vector contained in a circle Cr of radius 
R and area Ar around the origin. According to the classical definition the curl of the plastic distortion, the Burgers 
vector associated with k follows from Stokes theorem as b = nd^r which goes to a finite value as i? —>■ oo. If 

we instead evaluate the associated Burgers vector from the ’dislocation strain’ Cb, we get 

h = -[pb^D)/{A^p)C-^^^{MP)-J £ d^Kds,. (104) 

It is easily seen that the integral, rather than converging to a constant value, becomes exponentially small once the 
radius of the circle becomes much larger than £. It follows that, in such a theory, depending on the choice of the 
Burgers circuit an accumulation of excess dislocations may have no Burgers vector at all. We conclude that the back 
stress should better not be used to replace the long range dislocation interaction. 






